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Waveless ships in the low speed limit: 
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In the low-speed limit, a blunt ship modeled as two-dimensional semi- infinite body with a 
single corner can never he made waveless. This was the conclusion of the previous part of 



our work in Trinh et al. ( 2011 ), which focused on the Dagan & Tulin ( 1972 ) model of ship 



waves in the low speed limit. In this accompanying paper, we continue our investigations 
with the study of more general, piecewise-linear, or multi-cornered ships. The low-speed or 
low-Froude limit, coupled with techniques in exponential asymptotics allows us to derive 
explicit formulae relating the geometry of the hull to the form of the waves. Configurations 
with closely spaced corners present a non-trivial extension of the theory, and we present 
the general methodology for their study. Lastly, numerical computations of the nonlinear 
ship-wave problem are presented in order to confirm the analytical predictions. 
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1. Introduction 

The investigations in this paper are focused on the analysis of the low-speed, or low- 
Froud^ wave models proposed by Dagan & Tulin (1969 1972), in which blunt-bodied 
ships are studied in the context of potential flow and asymptotic expansions in powers 
of the Froude number. As a particularly interesting case that draws our attention, we 



recall the work of Farrow & Tuck (1995), who showed that by attaching a bulbous-like 



obstruction to an otherwise rectangular ship's stern, one could produce a dramatic effect 
on the production of transverse waves. As they reported in their paper: 

At this [Froude number], a rectangular stern generates waves with steepness 0.0855, 
whereas the stern with the downward-pointing bulb [. . . ] yields waves with steepness 
0.0119. It is clear that the addition of the downward-pointing bulb has had a dramatic 
effect on the downstream wave steepness, reducing it by a factor of 7.2, although it 
has still not eliminated the waves entirely. 

Our goal is to give an analytical criterion that explains why this phenomenon occurs; that 
is to say, what distinguishes the two ships, one with a bulb and one without a bulb, in the 
context of the 'slow-ship' approximation? The advantage of the slow-ship potential-flow 
approximation, is that allows us to directly relate the generation of waves to the shape of 
the ship's hull without the need for numerical simulations. 



f The Froude (draft) number represents the ratio between inertial and gravitational forces. 
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(a) 1-Hull 
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(b) 2-Hull 




Figure 1. Are any of these ships waveless? In all cases, the flow is from left to right and nodes 
indicate singularities in the analytic continuation. 



In addition to addressing the Farrow & Tuck ( 1995 ) issue, we are also interested in 



a more general question: in the low-speed limit, when a blunt ship is modeled as a 
two-dimensional semi-infinite body, can it ever be made waveless? These waveless or wave- 
minimisation questions in the context of the Dagan & Tulin ( 1972[) app roximation were 
studied by Vanden-Broeck & Tuck (1977), Vanden-Broeck et al. (1978), Madurasinghe 
& Tuck (1986), and Tuck ( 1991a|6 ) for ship hulls of varying geometries, and we are 
interested in continuing their line of inquiry. 

In the previous part of our work ( Trinh et aL||2011 ), we demonstrated that piecewise- 
linear hulls with a single, submerged corner can never be made wavelessif] For the case of 
piecewise-linear hulls with multiple corners, the answer to this question is not as clear. 
For example, are any of the eight hulls presented in Figure [l] waveless? If not, then 
which ones produce the smallest waves? For the case of potential flow over a submerged 
obstruction, waveless configurations are certainly possible, as was demonstrated by |Lustri] 
et al. (2012) and Hocking et al. (2012), but the same question for surface-piercing ships 



of general form remains open. Certainly, there are notable difficulties in studying this 



problem. For example, waveless ships were proposed by Tuck & Vanden-Broeck (1984) 



and Madurasinghe & Tuck (1986), but these were later refuted in the more comprehensive 



numerical study by Farrow & Tuck (1995), in which they showed that 



The free surface would at first sight appear to be waveless, hut on closer examination 
of the numerical data, there are very small waves present and they have a steepness 
of 1.5 X 10-3 



in reference to the bulbous hull in Tuck & Vanden-Broeck (1984). Our desire, then, is to 



study these issues in terms of the low-Froude asymptotic expansions. 

f Consequently, a free surface that attaches to a single-cornered bow at a stagnation point is 
not possible within the Dagan & Tuhn (19721) model. 



Exponential Asymptotics and Ship Waves 



Having progressed through the theory of Trinh et al. (2011), we know that at low 



Froude numbers, the waves generated by a ship become exponentially sniah and are thus 
invisible to a regular asymptotic expansion. The ineffectiveness of traditional asymptotics 



in capturing the low-speed limit was first remarked by Ogilvie ( 1968 



termed the Low-Speed Paradox. Techniques in exponential asymptotics ( Boyd||1998 ) allow 



1970 1 and later 



us to demonstrate the fact that these hidden waves are switched-on when the regular 
expansion is continued across critical curves {Stokes lines) in the complex plane; this 
process is known as the Stokes Phenomenon. Most important, however, is the valuable 
insight that these approximations give: an explicit formula that relates the shape of an 
arbitrary hull with its resultant waves. The question of wavelessness in the low-Froude 
limit is then simplified to examining whether the sum of the Stokes line contributions can 
ever be zero in regions far from the ship. 

The requisite background in exponential asymptotics can be found in [Trinh et al.\ 
(2011 ). The techniques we apply are based on the use of a factorial over power ansatz 



to capture the divergence of the asymptotic expansions, then optimal truncation and 
Stokes line smoothing to relate the late-order terms to the exponentially small waves [see 
for example, papers by Olde Daalhuis et al. (1995), Chapman et al. ( 1998[ ), and Trinh 



(2010a)]. Our paper also parallels the works of Chapman & Vanden-Broeck (2002 2006) 



and Trinh & Chapman (2013a|& ) on the application of exponential asymptotics to the 



study of gravity or capillary waves produced by flow over a submerged object. 



1.1. The role of low-Froude approximations and an outline of the paper 



It is important for us to mention that the low-Froude model of Dagan & Tulin ( 1972 ) is 



indeed a very idealised approximation for understanding the production of ship waves. 
Real stern and bow flows are very complex, and viscosity, turbulence, and necklace 
vorticies can all play an important role in the production of waves. We refer the reader 



to, for example, some of the numerical simulations of Grosenbaugh & Yeung ( 1989 ) and 



Yeung & Ananthakrishnan ( 1997 ) that demonstrate some of the complex dynamics that 
arise in ship flows once, for example, vorticity and viscosity are included. In ^of this 
paper, we shall return to discuss the caveats of the low-speed approximation. 

Ultimately, we are interested in obtaining analytical intuition about the connection 
between the ship's hull and the waves produced. The more usual routes towards analytical 
solutions assumes an asymptotically small geometry, which leads to the 'thin-ship', 'flat- 
ship' or 'streamline-ship' approximations; in such regimes, a waveless ship is impossible 
(see for example, Kotik fc Newman|1964 and Krein in Kostyukov|1968 ), but these theories 
say very little about the case of non-thin ships. Other examples of ship wave models can 
be studied, including the Kelvin-Neumann formulation in which the free surface condition 
is linearised about a steady uniform stream annd the boundary condition on the ship's 
hull is satisfled exactly, but generally these problems do require a degree of numerical 
computation. A discussion of such problems can be found in the book b y|Kuznetsov et al. 
(2002) and the review and discussion by Newman et al. (1991) (see e.g. Pagani & Pierotti 



2004 for more recent work on rigorous results applicable to non-slender geometries). We 
finally refer readers to the reviews by Tuck (1991a) and Tulin (2005) for a summary of 



the role played by low-Froude approximations, particularly in connection with problems 
in which we require asymptotic solutions that preserve the nonlinearity of the geometry. 
The paper will proceed as follows. First, the mathematical formulation of the ship-wave 
problem is briefly recapitulated in S|2J This is followed by the asymptotic analysis of the 
low-Froude problem in ^ which culminates with the derivation of explicit expressions for 
the wake of an arbitrary multi-cornered ship. From these analytical results, we explain 
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why certain classes of multi-cornered ships can never be made waveless in ^ then in f|5] 
these theoretical results are vindicated by comparisons with numerical computations. 



2. Mathematical formulation 

Consider steady, incompressible, irrotational, inviscid flow in the presence of gravity, past 
the semi-infinite body shown in Figure [2] which consists of a flat bottom and a piecewise 
linear front face. There is a uniform stream of speed C/ as x — >■ — c», and we assume that 
the flow attaches to sterr{f]at a stagnation point, x = and y = 0. 

The dimensional problem can be reposed in terms of a non-dimensional boundary- 
integral formulation in the potential (0, ^)-plane. The unknowns are the fluid speed 
q = q{(j>, V'), and streamline angles, 9 = 9{(p, -0), measured from the positive x-axis. The 
body and free-surface are given by the streamline ^ = 0, and we assume the free-surface 
(0 > 0) attaches to the hull (0 < 0) at = 0. The free-surface, with f/; = 0, is then 
obtained by solving a boundary-integral equation, coupled with Bernoulli's condition: 



\ogq= -4 



0{V) 



7I"7-oo V' - 



dip 



eq 



.dq 



(2.1a) 
(2.1&) 



where e = U^ / gK is related to the square of the Froude draft number with upstream flow 
U , and K is defined hy K = X]i=i ^i^ where N is the number of corners and 0* — ~Ki 
is the dimensional value of the potential at each of the corners. In this way, if = a^ for 
1 ^ i ^ N denotes the value of the potential at the corners in the non- dimensionalised 
problem, we have the property that X]i=i ^i ~ 1- The derivation of (2.1 1 can be found in 
Trinh et al. (2011 ), and the only difference is the choice of scaling for the Froude number. 

< the 



We shall refer to the A^-cornered piecewise-linear hull as an A^-huU. For 
geometry of the hull can be described by 



^^a 



J' 



(2.2) 



for afe < (/) < flfe+i where k — 1, . . . ,N , ajv+i = is the stagnation point, and TTCTfe is the 
exterior angle at the fcth corner (see Figure l2|. When iV = 2, we will sometimes also refer 
to the ship as a [(Ti,a2]-hull. 

In (2.1a|, we write the portion of the boundary integral over the negative real axis as 






'N+l 

n( 

.fe=i 



ak] 



log gs ((/-), 



(2.3) 



where ctm+i is defined according to the requirement that the free-surface approaches the 
uniform stream, with — > and g — > 1 as (/) — > oo; this gives 



C^TV+l 



N 



(2.4) 



The function q^ in (2.3) serves to distinguish the different sorts of piecewise-linear ships. 



f The reversible nature of potential flow implies that any stern flow can be reversed to bow 
flow, with the additional condition that there are no waves far upstream from the ship (the 
radiation condition). 
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Figure 2. Flow past a piecewise-linear A'^-huU. The A'^ corners of exterior angles Ttai, ■na2, ■ ■ ■ , t^ctn 
in the {x, y) plane (left) are mapped to u; = — ai, —02, . . . , — ajv in the complex potential plane 
(right). The stagnation point is u; = — ajv+i = 0. 



Note also that the product representation of the complex quantity qs can be ahernately 
derived by using a Schwartz-Christoffel mapping applied to the polygonal hull shape and 
a rigid, flat free surface. 



As explained in Trinh et al. (2011 1, in order to study the Stokes Phenomenon, we must 
'complexify' the free-surface, and thus allow g(0, 0) H' q{w) and 9{(j), 0) ^-> 9{w) to be 
complex functions of the complex potential, + iO i-> w. Analytically continuing (2.1a| 



and (2.1&) gives 



log q±ie = log qs{w)+ Jfe{w) 

2dq . . 

eq —— = — smt), 
aw 



(2.5a) 
(2.56) 



where the ± signs correspond to analytic continuation in the upper and lower-half (j) 
planes, respectively, and Jff denotes the Hilbert Transform operator on the free-surface. 



Jifeiw) 



if — w 



dip. 



3. Exponential asymptotics 

A single-cornered ship will always produce exponentially small waves in the low Froude 
limit, e — > 0; these waves are explained by the presence of a Stokes line which emerges 
from the singularity at the corner. For a multi-cornered ship, the analysis proceeds almost 



identically to Trinh et al. (2011), except now, each corner of the hull has the potential 



to produce Stokes lines and its own separate wave contribution. In this section, we shall 
briefly re-apply the methodology of the previous work, and provide the corresponding 
formulae for the case of an A^-hull. 



3.1. Late- order terms 

Here, we perform the asymptotic analysis which corresponds to analytically continuing the 
free-surface into the upper-half 0-plane; continuation into the lower-half plane produces a 
complex conjugate contribution, which we add to our results, later. 
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We begin by substituting the usual asymptotic expansions 



= J2 e"^" 



and 



E^>»' 



(3.1) 



n=0 



into (2.5a) and (2.5b) (with the + sign). In the hmit e — > 0, the leading-order solution is 
the rigid- wall flow of (2.31, 



^0 = 0, 
go = Qs 



N+l 



Y\_{w + ak) 



k=l 



while the 0{e) terms are 



h = -Qo 



dw 



qi = «9o-i \- qo-yiCOiiw). 



(3.2a) 
(3.2b) 

(3.3a) 
(3.3b) 



Notice that the leading-order solution, go in (3.2b), possesses a singularity at each of the 
corners, w — —a^. However, the solution at each subsequent order involves a derivative of 
the previous order, so we would thus expect that as n — >■ oo, the power of the singularity 
grows, and the asymptotic expansions (3.1) exhibit factorial over power divergence: 

^n - 2^ ^^+^^ and g„ - ^ ^^^^ , (3.4) 



k=l 



xt'^ 



k^l 



Xk 



where 7^ is complex constant, Q^ and Xk are functions of the complex potential w, and 
Xk^—oLk) = 0. For most of the analysis, however, we can simply choose one of the corners 
of interest and add the individual contributions at the end. 

The singularities, w = — a^, are located off the free surface, where the Hilbert Transform 
in (2.5a) is evaluated and so, as justified in Trinh et al. (2011), J^9n{w) is exponentially 
subdominant to the terms on the left-hand side for large n. At ©(e"), (2.5a) gives 



90 



iqiqn-i 



as n -> 00, and substitution into (2.56) gives the relevant terms at 0{e" 



qlq'n-i + iQn 



290909"-! + 29o9i9'i-2 - * 



9n-l9l 

90 



= 0. 



(3.5) 



(3.6) 



We substitute the ansatzes of ( |3.4[ ) into ( |3.6[ ), and this yields, at leading order as n — )• 00, 

(3.7) 



dx 
dw 



4' 



Using Xki—cik) = 0, we integrate this result, to give 



At the next order as 



Xk{w) -- 

> 00, we find that 

Qk{w) = ^-exp 
9o 



a, 9o('^) 



dip. 



3i 



91M 
9o(V5) 



dip 



(3.8) 



(3.9) 
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where A^ is constant, and w* is any point for which the integral is defined. FinaUy, (3.51 



allows us to relate Qk with 9^, 



using Gfc 

Afci 



4 



exp 



iQk/lo: so that 
91 M 



3i 



(laiv) 



dip 



With Xfc, Qfc, and 0^ now determined, we have thus derived the late-orders behaviour 



in (3.4), subject to the values of 7^ and A^; these must be determined by applying the 
method of matched asymptotics near the singularity, w = —a^- 

3.2. Stokes lines and the Stokes Phenomenon 



From Trinh et al. (2011), we know that the late-order terms (3.4) play a crucial role in 



determining the free-surface waves. Using the expression of Xk in (3.8 1, Stokes lines can be 



traced from each of the ship's corners, across which the Stokes Phenomenon necessitates 



the switching-on of waves. From Dingle (19731, these special lines are given by the points 
w; G C where 

3[XfcH]=0 and 3?[xfcH]^0. 



If we write go ~ Cfc(u; + a^) '^'° near w ~ —ak, then from (3.2b I, we have 



Cfc 



N+l 

i=i 



a-k) 



(3.10) 



and thus from (3.8), in the limit that w — > —a^, we have 



Xk 



4il + Sak) 



(W + Ok) 



l+3o-fc 



The condition that Xk{~o-k) = thus requires that ak > —1/3. In other words, for there 
to be a singularity, the local deviation of the corner must be greater than — 7r/3. This is a 
necessary (but not sufficient) condition for there to be a free-surface wave produced by 
the corner. In fact, a stronger condition for the existence of a Stokes line emerging on the 
relevant Riemann sheet can be derived. Since qe^^^ — u — iv, we can write 



arg(cfc) 



(3.11) 



for analytic continuation into the upper-half plane, where 9k G (— ti"; '"') is the angle of 
^, measured from the positive cc-axis (shown in Figure 2|. If we write 



the hull as w 
arg(w -I- Ofc) = 



i^k: then Stokes lines must leave at angles of 



J^fc 



SOk + 2m7r - 7r/2 
1 + 3afc 



(3.12) 



for ?7i G Z and we thus need v G (0, tt) in order for the line to emerge in the upper-half 
plane. The general requirements for a Stokes line to intersect the free-surface is a global 
function of the leading-order flow, but for most hulls, the requirement that v G (0, tt) with 
(3.12 ) is adequate. In any case, we will let J C {1, 2, . . . , iV -I- 1} denote those corners 



which have Stokes lines crossing the free-surface. 

As an example, consider Figure [3j which illustrates the Stokes lines for various A^- hulls, 
including a simple 2-hull, the 3- hull of Farrow & Tuck (1995), a bulbous 6-hull, and 



a step-like 9-hull. With the exception of a single configuration, the condition that a 
Stokes line emerges into the upper-half plane is enough to guarantee that it intersects 
the free-surface. The exception is with the 3-hull, for which the second singularity has a 
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9(w) 



Q{w) 



Q{w) 




I I I y I I i<^ I I >M< I I I V I I 1^1 I I w-r\/ i I I I I I I I I I I I I I I I I I I I I I I I I I I 
-1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 



5R(u;) 



Q{w) 



0.3-1 



0.2 



0.1 





I I I yT-WT-|M--r^7-yv-r-v.']--iM-Ty-rv-ry i i i | i i i { i i i { i i i { i i i | i i i | i i i 
1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 



5R(u;) 



Figure 3. From top to bottom: Stokes lines for the 2-hull, Farrow and Tuck's (19951 3-hull, 
the 9-hull, and the 6- hull shown before in Figure p] For the 2-hull and 6- hull, the corner angles 
diverge at ±7r/4; for the remaining hulls, the corner angles are all rectangular. 



Stokes line emerging at an angle of 1^2 = Btt/S, but which does not later encounter the 
free-surface. 

To derive the form of the exponentials that appear whenever a Stokes Line intersects 



the free-surface, we optimally truncate the asymptotic expansions (3.11, and examine the 
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remainder as the Stokes line is crossed. We let 



Af-l 

E 

n=0 



e"g„ + Sm, 



(3.13) 



with a similar expression for the series for 9. When J\f is chosen to be the optimal 
truncation point, the remainder 5V is found to be exponentially small, and by re-scaling 
near the Stokes line, it can be shown that a wave of the following form switches on: 

27ri 



g7fc 



Qfccxp 



Xfc 



(3.14) 



To (3.14), we must also include the complex conjugate due to the contributions from 
analytic continuation of the free-surface into the lower-half (/)-plane [see (2.5a)]. The sum 
of the two contributions is then 



Qcxp.fc 



47r 

( 

£7fc 



^<Qk exp 



e 



(3.15) 



with, of course, one such expression for every k €z J^. 

Thus, for any given arbitrary A^-hull with a geometry such that J7 is nonempty, the 
appearance of exponentially small waves is a necessary consequence of the divergent 
low-Froude problem; in order to check that such a ship can never be waveless, we need 
only verify that the sum of all the contributions incurred can never be zero, so that there 
is a non-zero wave amplitude far downstream. 

3.3. Wave formulae for N -hulls 



The constants 7^ and A^,, which appear in the final form of the waves (3.15) [the latter 
as a prefactor in Qk in (3.9)], can be determined by re-scaling w and q near each of the 
singularities, and then matching the leading-order nonlinear (inner) solutions with the 
late-order (outer) terms of (3.4). It can be shown (see (6.8) in |Trinh et aL|[20TT ) that 

6o-fc 

Ik 



l+3t7fe 



and 



where Ck is given by 



Kk 



2Ck (1 + Sa,.)"'' 



lim 



Ck = go(w*)exp f 3i / 

\ J w 



^00 r(?i + 7fc 



90 (v) 



dip 



The terms 



,^k are given by the recurrence relation, 
1, 

' ' 2ak 



oo.fc 



n-l 

E 

m=0 



l + 3afe 



^?n*rn— m— 1 



for n > 1. 



(3.16) 
(3.17) 

(3.18) 

(3.19) 
(3.20) 



We will often make reference to the limiting ratio in (|3.17|, so we define the function 

il,{(Tk) = hm 



L>n,k 



r(n 



Ik) 



(3.21) 



The value of 51(afe) only depends on the local divergence of the k corner, and its values 



are given in Trinh et al. (2011). Since fl ^ for all choices of the local angle Cfe, A^ is 
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also non-zero and this verifies that each of the \J'\ corners of an A^-hull must necessarily 
generate a non-zero wave on the free surface. 



With Qk given by ( |3.9[ ), A^ given by ( |3.17[ ) and ( |3.2ip , and arg(cfe) from (3.11), we 



have from (3.151 the result that 



Qcxp,k 



47r Icfcl^-^T*^ njak ) 
'e')''=2(l + 3CTfe)')''= q^ 



exp 



-5 3 



JfOi 



ak % 



exp 



J'/^Oi 



+ x + ^^ + (6-37fc)6ifc+5R 3 , 3 



(3.22) 



Then, for each k G J', we add the waves together, so that the total wave contribution 
after all the Stokes lines have been crossed is 



kej 



(3.23) 



4. The non-existence of waveless ships 

Let us see what would be needed to produce a waveless ship. 



The wave contributions (3.22) are written in terms of different denominators Xk (also 
referred to as the ^ singulants\ c.f. Dingle|[T973 p. 148). To make it easier to sum them we 
rewrite them in terms of the single singulant, xi- Note that 



Xiiw) = i / -^ +Xkiw), 

where in order for the integral to exist, we may have to avoid the intermediate corners by 
deforming the contour into the upper half plane. Consider now gcxp.fe in (3.22) when w 



is evaluated on the free-surface, that is, for w G M+. From the appendix of Trinh et al. 



(|2011 


), it was shown that 


Ufa 






exp 


\ „3 



ql{w)e. 



Moreover, the real part of Xi(w) comes from the pole at infinity, giving 



3ff(xi) = 3^ 



N 



(4.1) 



Putting these together in (|3.22|) we find, for w G 

Afe 



Q'cxp,/c 



<lli.w] 



exp 



3.^ 



cos 



3(xi(w)) 



3? 3 



Jf^i 



ai % 



where the dependence on k arises only through the constants 



Afc 



47re |cfc 



6-37fc 



eTfc 2(1 + 3crfe)T"= 



r2(o-fe)exp 



3 3 






exp 



*fc = -3ff 



ai % 



-3fi 3 






TT 

2 



TTTfe 

2 



** 



9o 



+ (6-37fc)0fc. 



4.1. On the two-cornered ship (2-hull) 



We have already shown in Trinh et al. (2011) that a single-cornered hull must produce 



waves. Therefore let us consider first the next simplest case of a 2-hull. For such a ship 
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to be waveless, both corners must generate Stokes lines which intersect the free surface, 
and the waves generated by each must exactly cancel; then, there will be a finite section 
of free surface containing waves, but no wavetrain at infinity (as happens in the case of 
capillary waves in^Chapman fc Vanden-Broeck|[2002). 



In order for the waves from the two corners to cancel, we need Ai = A2 . Now for a fixed 
value of e this may indeed be possible (and we give such an example in q5|, but what if 
we want the waves to vanish for all (small) values of e? Then, since each A^; is of the form 

Aie^^e-^^/% 

as e — > 0, we need the exponentials to be equal, the powers of e to be equal, and the 
prefactors to be equal. In order for the exponentials to be equal, we require 

'4) = 0. (4.2) 



ai 



% 



and since 



arg ( -3 ) = -36*^ for ak < w < Uk+i 



the only way (4.2 1 can hold is if 9i = 7r/3, so that q^ is real for —ai < w < —02. Now, for 
the algebraic factors of e to be equal, we require 71 =72, which implies cri = fT2. Thus 
the only possibility for a waveless 2-hull is for a ship with cti = (T2 = 1/3. To eliminate 
this final possibility we need to consider the prefactors. Since 

^ ^ \ 6i{w) w > 0, 
and Qq is real for —ai < w < —02. then 

The only remaining difference between the two prefactors is in Cfc. However, since 



c. 



,3 



(01-02) (01-02) 

the only way that we can have |ci| = |c2| is if oi = 02. Thus the two prefactors must be 
different, one corner always dominates the other one, and the waves can never cancel. 
Waveless ships with two corners are not possible. 

4.2. On general N -cornered ships (N-hulls) 

What can we say about more general ships? Let us take a general TV-hull with the following 
assumptions: suppose that all the Stokes lines intersect the free surface, that (ife > for 
each k, and that 0^ ^ 27r/3. In Figure 1, hulls (a) to (e) satisfy this requirement, whereas 
hulls (f) to (h) do not. Under these assumptions. Ok is monotonically increasing with fc, 
so that arg(l/5g) is monotonically decreasing. Thus the argument of the exponential 

1 ^ ( r"" d0 






is convex in fc, increasing while < 6k < 7r/3 and then decreasing while 7r/3 < Ok < 27r/3. 
Thus if Oj ^ 7r/3 for all j then we can see immediately that the waves generated at the 
corner fc such that 9k-i < 7r/3 < Ok exponentially dominate all the others. On the other 
hand, if Ok = 7r/3 then the waves from corners fc and fc + 1 have the same exponential 
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factor (as in the 2-hull case). If we further impose ak — <Jk+i then they have the same 
algebraic factor, and there is possibihty of wave canceUation if the prefactors are equaL 

Of course, even if we could get the prefactors to be equal, we still have to worry about 
the waves generated by all the other corners. In fact, even from these two corners there 
would be higher-order correction terms (both in the form ee~'^''^, e^e"'^''^, etc. and also in 
the form of a trans-series e~'^^'^, e~^'^/' etc.). Thus it does not seem to be worth pursuing 
the analysis further. However, even if we cannot get the waves to cancel exactly, we might 
expect a significant reduction in the amplitude of the waves in the case when leading-order 
cancellation occurs. 

This brings us to the natural question of whether the analysis we have presented may 
be used to design a hull to minimise the wave drag. Before we address this question, let 
us first demonstrate that the hulls shown in Figure 1 (f ) through (h) must generate waves 
on the free surface. 

We consider them in reverse order. The 9-hull shown in Figure 1(h) has \J'\ = 5 (as 
shown in Figure 3), with arg(l/gQ) alternating between zero and — 37r/2; thus the argument 
above can be used to show that the contribution from ai exponentially dominates all the 
others. The hull shown in Figure 1(g) has \J'\ = 4, with three positive angles cr and one 
negative angle. Thus the algebraic factor in the contribution from w — —a^ is different to 
the others, and those waves must always exist on the free surface. 

Finally let us consider the 3-hull shown in Figure 1(f), which is found in the work of 



Farrow Sz Tuck (1995), and for which the addition of a downward pointing bulb was 
shown to dramatically reduce the wave resistance compared to a rectangular ship. In this 
case, \J^\ = 1, and w = —03 is the only relevant corner, so there are always waves on the 
free surface. The principal effect of the bulb is to lower the usual singularity farther away 
from the free surface, thereby decreasing the amplitude of the waves. 

This last example not only highlights the difficulty in trying to minimise the wave drag, 
but also the advantage of our semi-analytic approach: we have gained considerable insight 
into the mechanism of wave production; from this, we can immediately see why Farrow & 
Tuck obtained the results that they did. 

In trying to design reduced-wave hulls, it is crucial to specify what exactly is the 
optimisation process. For example, if we simplify the hull of Farrow & Tuck to a 2-hull by 
reducing the width of the downward pointing bulb to zero, then we have one parameter, 
oi, over which we can optimise (since oi 4- 02 = 1). We find the smallest waves occur for 
Qi = 0.5, i.e. when there is no bulb. However, in our current nondimensionalisation, as 
we vary ai, both the depth of the hull and that of the bulb vary. If instead we fix the 
depth of the hull, and allow the depth of the bulb to vary, we find that the waves get 
smaller as the bulb gets deeper. On the other hand, if we fix the depth of the bulb, and 
allow the depth of the hull to vary, we again find that the smallest waves correspond to 
the hull depth being equal to the bulb depth, i.e. to there being no bulb. 



5. Numerical and asymptotic results for two-cornered hulls 



The numerical algorithms developed in Trinh et al. (2011) can be used to verify the 



asymptotic predictions. Here, we focus on the particular case of a 2-hull, which we refer 
to as a [(Ti, CT2]-hull; this is a ship with divergent corner-angles ci and CT2, and with 



leading-order flow given by (2.3), or 



^' {w + aiy^{w + a2y^' ^ ■ -* 



with ai + 02 = 1. 
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Figure 4. Solutions for the [0.5, 0.125]-hull (dashed hne) and [0.25, 0.25]-huU (solid line) at 
e = 2/3 and e = 1/3, respectively. Both ships have corners set at ai = 0.8 and a2 = 0.2. The 
solutions were computed using ALGORITHM A of Trinh et al. 12011 1 with n = 1000 and Acj> = 0.015 
for the former ship and n = 2000 and A0 = 0.015 for the latter. 



As we discussed in Trinh et al. (2011 ), the numerical computation of the stern problem 



at small values of e can be particularly difficult, and the culprit is the presence of the 
attachment singularity at w = 0, associated with a small boundary layer; this singularity 
is responsible for most of the numerical error. For hulls where the in-fluid attachment 
angle between the free-surface and body is less than 7r/3, a simple finite difference scheme 
based on the methods outlined in IVanden-Broeck & Tuck (19771 can be used, provided 



that we limit our search to waves larger than « 10~*. Figure |4j provides an example of 
solutions found using this method. 

The theory of ^ and fJ4] can be verified by comparing the analytical predictions with 
numerically computed wave amplitudes far from the ship. First, consider the effect of 
varying the Froude number on ships of a fixed geometry. This is shown in Figure [5] for 
three 2-hulls with their corners fixed with oi = 0.8 and a2 — 0.2. The individual cosine 
waves are calculated from (3.22) with qq ^ 1 downstream, and then the final amplitude 
computed using the sum (3.231 (see §4.5 in Trinh|20106 for additional details). The match 
between numerical and asymptotic solutions is quite good, and like the previous work, we 
remark that the results are applicable over a wide range of Froude numbers. 

Next, we would like to consider the effect of fixing the Froude number, but varying 
the positions of the corners. The problem, however, is that the interesting effects of this 
procedure are only seen at values of e much smaller than what we can achieve using the 
above numerical methods. In Appendix [Xj we present a slightly simplified version of the 
ship-wave problem (2.5a)-(2.5&) that preserves the asymptotic structure of the waves. 



but also enables us to compute numerical solutions to much higher accuracy. 

This simplified problem was used to create Figure |6] which shows the effect of varying the 
positions of the corners on the wave amplitude for a [j, j]-hull with e — 0.15, ai + a2 = 1, 
and values of 0.5 ^ oi ^ 1. The figure contains a number of interesting effects: the first 
is that the numerically computed wave amplitude shows a significant dip (an order of 
magnitude) near ai — 0.96, and that this effect is also captured fairly accurately in the 
asymptotic solution. 

The reason for the dip is that at this value of oi, the waves from the two corners exhibit 
partial destructive interference. However, from the of the previous section, we know that 
the waves generated by the corner at —02 should exponentially dominate those from — ai. 
How then are they cancelling each other? The answer is that the prefactor jci 



16-371 



IS over 



ten times larger than |c2| '^^] at this particular value of e, the difference in prefactors is 
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Figure 5. Numerical (dots) and asymptotic (solid) amplitudes of the downstream waves for a 
range of hull inclinations. In all cases, the corner points are fixed wit h ai — 0.8 and 02 = 0.2. The 
solutions were computed using ALGORITHM A of Trinh et al. (2011 1. The parameters used were 



the following: ® n = 1000, A0 = 0.04; © n = 1500, A(/> = O.OSTaiid ®n = 2000, Ac/) = 0.025 
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Figure 6. The numerical solution (solid) is plotted against the asymptotic approximation 
(dashed) for the simplified nonlinear problem of ijAl The ship is a [|, j]-hull with ai + 02 = 1 
and e = 0.15. The dotted line indicates the one-cornered approximation for a rectangular hull. 



enough to compensate for the difference in exponentials, since 



^ 



1I 



is not very large. Thus, the value of oi at which cancellation occurs depends on e; for 
somewhat smaller of e the Stokes line from w = — 02 does indeed dominate and no 
cancellation is possible. This indicates that it should be possible to design hulls with 
reduced wave drag at a particular speed (Froude number) . It is also reassuring to note that 
including the leading-order term from each of our exponentially-small waves captures the 
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behaviour of the solution very weU, even though formally, one of the terms is exponentially 
subdominant. 

The second effect illustrated by Figure [6] is the divergence between our asymptotic 
expansion and the numerical solution for values of ai close to 0.5. The reason for the 
divergence is that the prefactors Ci and C2 are singular as the corners approach each other; 
our analysis in Sj3] in fact, implicitly assumes that the corners of the ship are spaced 
sufficiently far from one another. 

To be more specific, in order to determine the constants "fk and A^ in (3.161 and (3.171 



in the previous analysis, the outer solution was required to match a nonlinear inner 
solution. The size of this inner region can be derived by observing where the breakdown 
in the outer expansion (3.11 first occurs, i.e. where eqi ~ qo. From (3.2bl and (3.3b), the 



required re-scaling of w near a singularity at — a^ can be seen to be 



Qk = O (ei+3" 



(5.2) 



Thus, if two (or more) singularities are spaced within a distance of (5.2) apart, then the 
previous asymptotic methodology breaks down as e — > 0. If we combine the two corners 
into a single corner of angle 7r/2, we find the wave amplitude is given by the dotted line 
in Figure [6J This approximation clearly works well for ai close to 0.5. 

A uniform approximation would need to smoothly match with the one-cornered ap- 
proximation at one end, and the (separated) two-cornered analysis at the other, and thus 
bridges the dashed and dotted lines in Figure |6) Such an approximation requires us to 
consider the distinguished limit in which the corners of the ship are allowed to approach 
one another as e — >■ 0. This is similar to the situation in ( Chapman fc Vanden-Broeck|2006 
Appendix B) where the asymptotic solutions for flow over a rectangular step in a channel 
was analysed in the case where both corners of the step lie in the same 'inner region'. For 
the case of the multi-cornered ship, the details of this analysis are very technical, and 
will be published in a future paper. 



6. Discussion 

In the end, what is the definitive answer to the conundrum of existence and non-existence 



of waveless ships? Unlike our results for the single-cornered ships of Trinh et al. (2011 ) 



there does not seem to be a simple answer to this question, applicable to the most general 
piecewise-linear hulls. Despite this, however, we have offered several new insights into 
the study of ship-wave resistance: we have offered explicit formulae for the computation 
of waves given the shape of the ship's hull; we have offered simple interpretations of the 
production of such waves in terms of Stokes line crossings and the Stokes Phenomenon; 
and perhaps most importantly, we offered a methodology which, given specific ships, 
provides an immediate and intuitive understanding of the effect of the body on the 
free-surface. 

In the previous work, we highlighted the importance of distinguishing between local 
and global properties of the analysis. Consider the factorial over power divergence of 



the asymptotic series in (3.4), or the emergence-conditions of Stokes lines in (3.12), or 



the numerically-determined pre-factor, 51^, in (3.21) — these are all local properties of 



the problem; indeed their derivation only depends on the behaviour of the asymptotic 
solution near the relevant singularities. These local properties, we understand well. 

In contrast, many global questions remain unanswered. For example, given a ship, 
represented by qa, what are the necessary and sufficient conditions for Stokes lines to 
intersect the free surface? Or perhaps more difhcult: what are the necessary restrictions 
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Figure 7. Stokes lines for the smoothed huh in ( |6.1[ ) with ai = 0.8 and a2 = 0.5. The Stokes 
hne leaves tangentially along the boundary, but later intersects the free-surface. 



on qq so that total phase cancellation occurs? We have provided a few preliminary results 
on this global problem, but a more exhaustive analysis of these issues remains an open 
problem. 

Naturally, our study of ship waves would be incomplete without a theory applicable for 
smooth hulls, with the eventual goal of addressing the well-known technique of using a 
bulb to reduce the wave resistance of a ship ( Baba|[l976 ) . However, the difficulty here is 
that analytic continuation is an ill-posed process and small perturbations in the shape of 
a hull can have large effects on the associated singularities — a unified theory for arbitrary 
ship geometries will likely prove difficult, if not downright impossible. 

Perhaps, then, we should only consider specific classes of smooth ships. Ship waves 
associated with continuous geometries have been considered in the numerical work of 



Tuck & Vanden-Broeck (1984), Madurasinghe (1988), and Farrow & Tuckl ( 1995 ), where 



there, the hulls are specified by piecewise-entire functions. For example, Farrow &: Tuck| 
( 1995 ) consider the family of hulls given by 





A{w - 



l)iw + b) 



■J! (tU+1) 

2 (1-fc) 



for w e (— oo, —1) 
for w E (—1, —6) 
for w e (-6, 0) 



which, given parameters A and 6, provides a ship consisting of a horizontal bottom and 
a vertical line, joined by a rounded section; ^ > yields rounded corners and A < 
yields bulbous sterns. The key is that if we restrict ourselves to classes of ships given by 
piecewise entire functions, then the complex singularities must be located at the points 
joining each piece. As a simpler example, we may consider the ship with 




for w G (— cx),— fli) 
for w e (— fli,— a2) 
for w e (—02,0), 



(6.1) 



which is similar to the vertically-faced one-cornered ships studied previously, but with 
a rounded edge. Analysis of the Stokes lines shows that the relevant line emerges from 
w = —a2', this is shown in Figure [7] The study of these piecewise-entire ships is the 
subject of ongoing investigation. 

A similar direction for research is towards the development of a low-Froude asymptotic 
theory for flows past three-dimensional, full-bodied ships. This builds upon the works 



of, for example, Keller (19791 and 'Brandsma & Hermans (1985), who apply geometrical 



ray theory to the case of streamline (thin) ships. In theory, the interpretation we have 
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presented in this paper of free-surface waves arising due to Stokes line crossings is still 
valid in three dimensions, except now, singularities are associated with Stokes surfaces 
rather than lines. In practice, however, the analysis is complicated due to the loss of 
complex variable techniques. We refer the reader to the work of Chapman & Mortimer 



( 2005 1, which provides a first step towards extensions of exponential asymptotics to partial 



differential equations. 

In addition to our study, which solely focuses on the low-Froude model of [Pagan fc 
Tulin (19721, it is important for us to question the place of these simplified mathematical 



models in terms of the bigger picture: that which includes the effects of vorticity, viscosity, 
and time-dependence in ship-wave interactions. As we elucidated in the introduction, 
numerical work (as particular examples, see Grosenbaugh fc Yeung|[l989 and Yeung fc 
Ananthakrishnan]|1997 1 show that in practice, these neglected effects can have significant 
roles in the production of waves. Extended discussions of the role of low-Froude theories 
appear in [ Tuck fc Vanden-Broeck] ( p84l p. 301), |TiIck|fl991a[ ), and |Mm| (|2005j) . Thus, 



while it is certainly true that in order to obtain analytical approximations directly relating 
ship geometries to free-surface waves, the low-Froude approximation provides enormous 
simplification, we hope that similar analytical theories can be developed which include a 
more complete host of effects. 
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Appendix A. The simplified non-linear problem 

The full problem in (2.5a)-(f2.5&) can be studied using the methods we develop here, 



but can also work with a simpler problem that nevertheless contains all of the the key 
components. The reason for this simplification is that, in order to verify the asymptotic 
analysis in the regime where the ship's corners are closely spaced, wave amplitudes must 
be computed to five or six digits of precision — otherwise, the fine effects of adjusting 
the ship's geometry are easily missed; this precision can only be easily achieved for the 
simpler problem, which we now derive. 
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As we know, when exponentially small terms are sought from (2.5a), the integral term, 



Jif9, only serves a minor role throughout the analysis. If we return to the derivation of the 



late-orders ansatz (3.4), we recall that the subdominance of J^6 as n — )■ oo ensures that 



it plays no part in the derivation of Xk- In fact, the only role of the Hilbert Transform is 



to change the expression for qi in (3.3b). Consequently, in the final form of the waves 



Qoxp.fe in (3.22), the presence of J^9i only serves to change the amplitude coefficient and 
the phase shift by an 0(1) amount. 

Therefore, the salient features of the problem can still be retained if we use log q±i6 = 



logQs instead of (2.5a); this way, we simplify the full problem in (2.5a) to (2.56) to a 



single nonlinear differential equation in q. Analytic continuation into the upper-half plane, 
and substituting i9 



\og{qs/q) into (2.5&) gives 
„3dq , 



e9s9 



dw 



Is 



= 0, 



which can be solved subject to the single boundary condition q{0) 



convenient to work under the substitution u(w) 

du 



eqsU 



dw 



q (w), where we have 



0, 



(Al) 
0. It is more 

(A 2) 



as a simplified nonlinear model of ship waves. Simplifications of the boundary integral 
problem (2.5a) and (2.56) were also proposed in Tuck ( 1991a|6|, b ut there, the simplifica- 
tions were argued based on behavioural requirements. Here, ( A 1 ) is a justified reduction 
based on the e — >■ limit. 

Notice that in this new problem, we chose to analytically continue into the upper 



half-w-plane, and thus the exponentially small waves of ( A 1 ) will possess both a real and 



imaginary part. If we wish, we can mirror the analysis for the lower half-w-plane and add 



the complex conjugate as we did before for (3.15) 



However, it is somewhat simpler to examine ( A 1 ) as a problem on its own; thus we shall 



only concern ourselves with studying the real component of the solution to ( A 1 ) , which 
we write q^^p — 5ft(<?oxp)- Now instead of (3.22), the form of the waves for the simplified 
problem (with well-separated corners) is given by 



ycxp.fc 



27r \ck 



|6-37fc 



ri(crfc) 



e^>'2{l + 3ak) 



Ik 



ll 



exp 



3ff(x* 



cos 



3(Xfc) 



TTTfc 

2 



(6 - 3jk)0k 



(A3) 



which is effectively (3.22 ) with Jf = and without a phase shift of tt/2. The reduction 



by a factor of 2 in (A3) compared to (3.22) occurs because there is no need to add the 



complex conjugate wave contribution. Analytical and numerical results for the simplified 



et al. 



nonlinear pr oblem of (|A 1|) in the context of a one-cornered ship can be found in 

1 r 

4 



(2011 ), whereas we have already discussed the numerical solution of the 



Trinh 



-hull 



for the simplified problem in ^ and Figure [6j 



